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(57) Abstract: A method and a 
system for controlling tlic wheel 
slip (A.) on a vehicle is proposed, 
comprising the steps of obtaining at 
least one dynamical state parameter 
of the vehicle, determining (G(s)) 
a slip reference value indicative 
of the wheel slip at a determined 
driveline torque on the basis of a 
dynamic model of the relevant vehicle 
components, and controUing (R(s)) 
the driveline torque (Ts) of the vehicle 
by determining a torque manipulation 
signal based on the dynamical state 
parameters and the slip reference 



value. 
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METSOD, SYSTEM AND COMPUTER PROGRAM PRODUCT 
FOR WHEEL SUP CONTROL 

FIELD OF THE INVENTION 

The present invention relates to wheel speed control of a vehicle, and more particularly to a 
method, a system and a computer program product for carrying out a spin or slip control on 
one or more vehicle wheels or other type of driven wheels. 

BACKGROUND OF THE INVENTION 

The automobile industry is one of the biggest and most competitive industries today. There are 
many actors and an abundance of different cars aimed at a highly diversified market and great 
efforts are made to improve performance and safety in the products without reducing profits. 
One strategy for doing this is to develop systems that easily can be implemented in different 
cars and adapted to the customers' wishes. An example of this is the so calle d ABS-s ystems 
(anti-blocking system) that often are developed by one specific company and then sold to many 
different car manufiictures but properly adjusted to differences in vehicle characteristics such as 
vehicle mass and engine performance and customer preference such as desired comfort. 
Therefore, there is a need to develop a Vehicle Stability System with this strategy in mmd. 

A vehicle is e?q>osed to many different forces during driving. The sum of those forces will 
dedde how the vehicle will behave. The driver controls the vehicle trajectory by controlUng the 
engjne torque through the accelerator, changing the steering angle, braking and/or using the 
clutch. However, the effects of the actions taken by the driver will depend on environmental 
conditions such as tire-ground traction, slopes and temperature. This unplies that the driver has 
to make a judgement of the environment before taking action. Sometimes this judgement is 
difficult to make correctly and this is where the vehicle stability systems comes in play. The 
most familiar Vehicle Stability System is probably t he ABS-brakes which is used for o ptimisii^g 
the braking force depending on the groimd conditions independently (the ABS-system cannot 
increase the braking force,' only reduce it) of the braking force requested by the driver. 
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Other vehicle properties that are of interest to control are slip, slide and skid. Slip is essentially 
defined as the speed difference between the vehicle and the driven whed and its value will 
decide the longitudinal and lateral forces fi-om the wheels acting on the vehicle. Slide is defined 
as the vehicles behaviour vM& cornering {under- and oversteer) and is partly dependent of the 
dip. Skid is here used to denote the unstable state when all the vehicle tires has lost their 
traction ability and the vehicle is out of control. There are many different systems available 
today on the market that tries to control those properties. Examples pertaining to the prior art 
are the Volvo systems TRACS, DSA and DSTC. They are quite different in their configuration 
and will not be finrther discussed here. 

One purpose of vehicle stability ffystems such as an anti-slip system is to control the wheel 
speeds of a vehide and by doing so, achieve optimal accderation and handling. This is done in 
the prior art systems with enq)irically devdoped controllers, which are usually implemented as 
large finite state madiines based on look-up tables usmg the engme and/or the brakes as 
actuator. Since the parameters of the anti-slip systems normally are configured for different 
driving conditions and vehicle specific properties, the implementation of the systems in new 
platforms has shown to be time-consuming and costly. 

OBJECT OF TBDE INVENTION 

It is an object of the present invention to provide an improved method, computer program 
product for carrying out the method and a system for whed speed control in vehicles. It is a 
fiirther object to provide a method and a system for tire slip control with the aim to increase 
flexibility and deo^ease costs and time for inq)lementation in different types of vehicles. 

SUMMARY OF TBDE INVENTION 

The invention provides a method for controlling the wheel slip on a vehicle, in particular 
carrying out anti-slip, anti-spin or ABS control, comprising the steps of obtaining at least one 
dynamical state parameter of the vehicle, determining a sUp reference value indicative of ti^ 
wheel slip at a determined driveline torque on the basis of a dynamic modd of the rdevant 
vdiicle components, and controlling the driveline torque of the vehicle by detenninii^ a torque 
manipulation signal based on the dynamical state parameters and the slip reference value. 



wo 2004/067307 



- 3- 



PCT/£P2003/000807 



According to another aspect, the invention provides a method for determinmg the traction 
force of a tire on a vehicle, in particular in combination with one of the preceding methods, 
based on a slip-fiiction model according to the formula = Mi^Wz / where is a fiiction 

value, A is a tire slip value, Ftrac is a traction force, i^,is the vertical load on a wheel and /j(X) is 
a slip curve according to the slip-fiiction model. 

According to still another aspect, the invention provides a system for controlling the wheel slip 
on a vehicle, comprising a device for obtaining at least one dynamical state parameter of the 
vehicle, a device for determining a slip reference value indicative of the wheel slip at a 
determined driveline torque on the basis of a dynantiic model of the relevant vehicle 
componmts, and a device for controlling the driveline torque of the vehicle by determining a 
torque manipulation signal based on the dynamical state parameters and the slip reference 
value. 

According to yet another aspect, the invention provides a computer program product including 
program code for carrying out a digital signal processing method, when executed on a 
compute system, for controlling the wheel slip on a vehicle, comprising a device for obtaining 
at least one dynamical state parameter of the vehicle, a device for determining a slip reference 
value indicative of the wheel slip at a determmed driveline torque on the basis of a dynamic 
model of the relevant vehicle connponents, and a device for controlling the driveline torque of 
the vehicle by detenmning a torque manipulation signal based on the dynamical state 
parameters and the slip reference vahie. 

The invention is applicable for and adaptable to all kinds of vehicles carried on pneumatic tires, 
such as cars, motorcycles, lorries or trucks having any number of wheels in accordance with 
the different embodiments and fiirther aspects described in the detailed description. 

BRIEF DESCRIPTEON OF THE DRAWINGS 

The present invention will be further explained by means of exemplifying embodiments in 
conjimction with the accompanying drawings, in which the figures schematically show: 
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FIG. 1 Diagram of a complete system enviromnent 

FIG. 2 Diagram iOustrotiiig the internal signal routing in vehicle. (The dash-dotted 
arrows indicate fhe route while ui»ng xPC Target). 

FIG. 3 Table of torque reduction actions, which are performed in a fixel blocking 

component. 

FIG. 4 Table compriiang principal patterns for the different actions in a fuel blocking 

conqponent. Non-blocked cylinders are indicated by (0), blocked cylinders are 
indicated by (1). 

FIG. S Table showing an example of an action sequence of a fuel blocking component. 

Performed sequence (bold) and planned sequence (r^;ular). The pattern for the 
applied action is underlined. 

FIG. 6 Schematic view of a spark ignition engine. 

FIG. 7 Diagram of the engme delay as function of the engine speed. 

FIG. 8 Schematic view of a three state model of a vehicle driveline. 

FIG. 9 Schematic view of a tire seen from above with slip-angle a. 

FIG. 10 Diagram of the slip curve with linear (a), non-linear (b) and unstable (c) r^ons 
(thermal (solid) and isothermal (dashed)). The traction parameter □ is 
diagramted as fimction of the slip parameter □. 

FIG. 1 1 Schematic description of a single wheel model. 
FIG. 12 Schematic description of a stiff driveline model. 
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FIG. 13 Diagram of the relationship between the real slip curve (dashdotted) and the 
approximation used in the controller (solid). 

FIG. 14 Validation of Gnom(s), Model output (dash-dotted) compared with measured 
data of slip (solid). 

FIG. 15 Validation of Gatsfs), Modd output (dash-dotted) compared with measure data 
of 2fv (solid). 



FIG. 1 6 Block diagram of the principal structure of a controller 



FIG. 17 Diagram of the step response for a closed loop system controlled with a IMC 
designed controller. Low □ (dash-dotted), "good" □ (solid), high □ (dashed). 

FIG. 18 Schematic diagram of the Ackerman compensation. 



HG. 19 Diagrams depicting the system identification. Measured output from vehicle 
(solid), model output (dash -dotted). The vehicle is hanging in mid-air to 
eliminate groimd influence. The values on axes are raw binary data. 



FIG. 20 Diagram of the sUp curve for gravel. 

FIG. 21 Diagram showing an estimation of kiip while driving fi-om asphalt onto gravel 
with winter tires. 



FIG. 22 Diagram showing the slip as fimction of time for a fiiU throttle take-off with 
DSA, initial engine speed 3000 RPM. 

FIG. 23 Diagram showmg the slip as function of time for a full throttle take-off with 
DSA, initial engme speed 3000 RPM. 

FIG. 24 Diagram showing the slip as function of time for a foil throttie take-off with 
NSC 2. Initial engine speed 3000 RPM. 
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FIG. 25 Diagram showing the slip as function of time for a full throttle take-off vnth 
NSC 2. Initial engine speed 3000 RPM. 

FIG. 26 Diagram showing the slip as function of time for a fiill throttle acceleration with 
DSA from 30 to 60 km/h. 

FIG. 27 Diagram showing the slip as function of time for a fiill throttle acceleration with 
DSA from 30 to 60 km/h. 

FIG. 28 Diagram showing the slip as function of time for a fiill throttle acceleration with 
NSC 2 from 3 0 to 60 km/h. 

FIG. 29 Diagram showing the slip as fimction of time for a full throttle acceleration with 
NSC 2 from 30 to 60 km/h. 

FIG. 30 Diagram showmg the slip as fimction of time for cornering with DSA. 

FIG. 3 1 Diagram showing the slip as function of time for cornering with DSA 

FIG. 32 Diagram showing the slip as function of time for cornering with NSC 2. 

FIG. 33 Diagram showing the slip as fimcdon of time for cornering with NSC 2. 

FIG. 34 Diagram showing the slip as function of time in an acceleration lane with DSA 

FIG. 35 Diagram showing the slip as fimction of time in an acceleration lane with the 
controller of the invention. 

FIG. 36 Block diagram of a further embodiment of the control structure of the invention. 

FIG. 37 Diagram of a fiiction ellipse showing how the total fidction force can be split up 
between lateral and longitudinal forces. 
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Diagram showing how the slip reference s^^ is computed from 0^. 

Three diagrams, where the first diagram shows the noise-free and noisy slip in 
an anti-slip state, the second diagram shows a true slip curve (solid), an initial 
guess (dotted) and a final estimated curve (dashed) using noisy slip 
measurements, and the third diagram shows the parameter convergence. 

Four diagrams, where the first diagram shows the convergence to iW^/wfrom 
different initial states, the second diagram shows the convergence to r,n from 
different initial states, the third diagram shows initial slip curves, and the fi>urth 
diagram shows final estimated slip curves. 

The left diagram shows the trajectoty fi^r s and ju on left and right driven wheel, 
respectively, during a drive with fiiUy developed spin. The first subdiagram fi^r 
static data, and the second subdiagram for LP-filtered //. The right diagrams 
show time diagrams of s and /i and the estimated slip curve. Three different test 
drives are compared in the different rows. 

DETAILED DESCRIPIION OF TOE EMBODIMENTS 

Many of the advantages of the present invention will be readily appreciated and become better 
understood by reference to the following detailed description when considered in connection 
with the accompanying drawings. However, before proceeding fiirther with the detailed 
description a few general obs^ations are made to improve imderstanding of the present 
invration. 

[Model-based Control] 

The known anti-slip systems on the market today are typically empirically developed 
controllers based on look-up tables with large numbers of parameters. The parameters of the 
controllers are normally dependent of the operating conditions such as gear, fiiction and 



nG.38 
FIG. 39 

FIG. 40 

FIG. 41 
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velodty. Whm sudi a system is to be ad^ted to a new platform an extensive testing 
procedure is required to tune the controller with regard to the new vehicles properties. 

A particular advantage of using model-based control is the flexibility of the resulting control 
system. The idea is to develop a general parametric control system that can run on any 
comparable platform of known characteristics without any significant loss of performance. In 
model-based control the number of parameters to be trimmed is substantially lower and thus it 
is easier and less time consuming to adjust the controller to a new platform. 

[Strategy] 

When developing model-based controllers it is necessary to first make a survey of the 
behaviour of the system to be controlled. This includes not only its physical and dynamic 
properties but also how the system environment is configured in a broader sense. The 
questions to be aiiswered are among others: 

• What types of sensors are available? 

• How will the data be transferred? 

• What is the limitation in processor power for the controller? 

• What type of actuator strategy should be used? 

• What other systems is the controller dependent of? 

• In what environment is the system meant to fimction? 

• What are the cases to handle? 

When these questions have been answered a scheme for the complete system environment can 
be established as illustrated in FIG. 1. The anti-slip system can be divided in two subsystems, 
the controller and the reference value calculator. The suggested controller was developed using 
Intemal Model Control (IMC). The reference value calculator was of secondary interest in this 
study and the version used in the implementation was just for evaluating the p^ormance of 
the controll^. 

When the complete system environment has been established the components can be 
characterised and moddled in a simulation environment to reproduce the behaviour of the real 
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system. From this mathematical model a controller can be derived by using a combination of 
model-based techniques and logic as described hereinafter. This controller can be evaluated in 
a simulated environment and if it seems to have promising qualities then tested on a real 
system. 

When the controller is unplemented it is often necessary to change some of its parameters to 
compensate for un-modeUed dynamics and other disturbances, especially if the controlled 
system is hard to describe. This implies that it is often rewarding to begin with quite simple 
methods where changes in parameters and configuration are easily performed and interpreted. 

[The Volvo ¥40 T4 Environment] 

In this embodiment a Volvo V40 T4 was used as platform during development of the 
controller and the system environment was configured as in FIG. 2. This embodiment is an 
approach to implement model-based anti-slip control and is mainly limited by the scarce 
amount of available signals in the V40. Thus, the system is to be kept simple to make it 
possible to evaluate its performance under real driving conditions. In the future, as the amount 
of available signal in cars rise through the use of CAN-buses and additional sensors such as 
accelerometers, the anti-slip system could be upgraded to take these in count. 

A Model-based Anti-slip Controller (MASC) according to the present embodiment has the 
potential to reduce complexity and implementation time, and makeUse of possible information 
fi^om other systems (such as fiiction estimation) and non-standard sensors. The idea is to use a 
model-based control system, which in this case include: a dynamic model of the tires, driveline 
and car body, a state estimator (ECalman filter) to compute the vehicle states and a cascade 
coupled controller to reduce the engine torque during fully developed spin. 

Here the actuator has been rqplaced by the actuator control system, the Engine Electronic 
Control Unit (EECU), since it seems more logical to implement the engine as a component in 
the vehicle. The term NSC is here used to denote the ECU for the anti-sUp system. In the 
present embodiment of the system the anti-slip algorithm is be placed here. However, while 
developing the controller the NSC only served as a communication node for all data traffic 
between the vehicle and a rapid-proto^ing tool, xPC Target, running on an external PC. The 
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EECU, NSC, sensors and some actuator related subjects are described hereinafter. The main 
part of the system environment is the actual vehicle or G(s) and its operational environment. 

[ENGINE ECU AND SENSORS] 

5 

This section deals with the parts of the system environment that concerns the generation of the 
torque reduction requested by the controller and information about the various electronic 
devices in the vehicle used for data collecting, transfer and treatment! 

10 [Engine ECU -EECU] 

The EECU receives the control signal from the anti-slip system and its task is to force the 
engine to reduce the torque. There are three basic methods for doing this, by controlling the 
throttle-angle, by controlling the ignition angje or by blocking the fuel to the cylinders. In the 

15 presented embodiments, fiiel blocking is preferably used for controlling the driveline torque 
e.g. engine torque but the controller could use any (or combinations) of the suggested 
methods, although its performance might be reduced. Another powerfiil method of controlling 
the slip woiild be to use the brakes to bring down the wheel speeds. However, this is quite 
hazardous and difiBcult to implement and it will also lead to that the anti-sUp system would 

20 become classified as a safety system with much higher specification requirements to meet. 

[Fuel Blocking] 

The fiiel blocking control technique is based on series of so-called actions. An action is named 
25 after how many cylinders that are blocked in an unagmed sequence of cylinders (16 qdinders 
or eight en^e revolutions in the case of the V40). This implies that this technique does not 
allow the engine to deliver the desbred torque reduction in a continuous form. The torque 
reduction has to be quantified to correspond to a specific action at the current engme torque 
requested by the driver, r^wvcr. A table, which lists the torque reduction percentage for each 
30 action is given in HG. 3. The procedure is easiest e)cplained by a simple example. 

Example: Say that Tdriver equals 200 Nm and that the controller suggests a torque 
reduction of 83 Nm. This means that the desired torque output &om the engine 
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Te equals 1 17 Nm, that is 58.5 % of Tdri^, The closest vahie in the action table 
(FIG. 3) is action 3 that gives 60.1 % so action 3 is sent to the en^e control 
system. 

The values found in the table of FIG. 3 were empirically estimated by measuring the vehicle 
acceleration at fiill throttie with diflFerent fixed actions set. The values for action 14-16 are very 
uncertain and depend on the prevailing throttie position (a small throttie angle probably gives a 
negative torque due to low pressure in the intake). 

ft 

Each action is associated with a certain pattern of blocked cylinders. When a sequence of 
actions is sent to the EECU their patterns are careMy fitted to each other to keep the same 
cylinders blocked as long as possble. This technique is referred to as cascade control. It is well 
known that unnecessary toggling between blocking and non-blocking for a cylinder will lead to 
a sub-optimal air-fiiel nux (stochiometric ratio) with undesired phenomena such as veiy hot 
exhaust gases and wear on the catalyst as consequence. The principal patterns for the different 
action are presented in the table of FIG. 4. 

If a sequence of actions is sent to the EECU the pattern-fitting algorithm results in a somewhat 
strange sequence of blocked and non-blocked cylinders as shown in the table of FIG. 5. The 
resulting pattern implies that the engine only can produce the desired torque in a broad 
statistical sense. At low engine ^eeds it will be possible to send actions to the EECU even 
&ster than th^ are possible to execute due to the large engine delay later described. This will 
lead to even more disturbances in the blocking sequence and more imcertainty regarding the 
resulting engine torque. 

[Sensors & Signals] 

The sensors for measuring dynamical state parameters of the vehicle such as wheel speed, 
engine torque, vehicle speed, vehicle traction, etc., are all ready existing sensors in the standard 
edition of the Volvo V40. The signals used are found in ahnost any vehicle on the market 
today. 



[ABS] 
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The wheel speeds were measured by the ABS-system. The sampling was perfonned at a rate of 
200 Hz and the signal could be seen as quite accurate. However, the signal will be rather noisy 
due to the wheels exposure to environmental changes. 

[EECU] 

The EECU measures the engine torque and speed and combines them in a Pulse Width and 
Frequency Modulated-signal sent at 50 Hz. The torque signal is the most uncertain due to the 
fact that its value is taken from a static map of steady-state torques. When the driver demands 
rapid changes in torque the signal will be &r from accurate. In addition to this the signal will be 
multiplied with the values in table of HG. 3 when there are actions present. 

[Filtering] 

Before the signals are used for calculations they have to be filtered in order to get rid of 
unwanted noise from the measuring process. It is assumed that the major part of this noise is of 
higher frequencies, at least it is safe to say that the "real" signal do not contain very high 
frequencies (this is normally the case with physical systems). Therefore a low-pass filter is 
apphed to reduce the noise influence. When filtering signals there is always a trade-ofT to be 
made between the desire to have a short rise time for the system in order to detect rapid 
changes (here e.g. ice patches), and the necessity to block out high frequency disturbances. In 
order to achieve this an Infinite Lnpulse Response (HR) filter is the natural choice. The down 
sdde with the IIR filters is that they all have a more or less non-linear phase. The Non-linear 
phase tends to delay signals differently as they pass through the filter depending on thdr 
firequendes. This phenomenon can give rise to strange impulse responses (for fijrther 
infi>nnation see [10] or any signal processing textbook) and is thus not deinrable. In this 
embodiment of the invention a Butterworth filter has been chosen because of its not so non- 
linear phase. 

[Data Treatment] 



[NSC] 
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For nonnal use all the calculations are perfonned in the NSC. However, during the 
development of the controller it was only used as a conununication node for the information 
exchange between the ABS, EECU and the xPC Target system. 

[xPC Target] 

xPC Target is a toolbox developed by MathWorks for rapid prototyping of real time 
applications. It gives the possibility to have SIMULENK® models on a PC (Host) and from this 
model generate C-code. The code could then be compiled and downloaded into another PC 
(Target) running a simple real-time kernel. This gives the possibility to make changes in system 
configuration and parameter settings and easily observe how various signals flow through the 
model. The SIMUUNK® model is ea^y modified and downloaded to the target computer and 
thus several different strategies can quickly be tested in the real system. 

In order for this to woric property, the communication between the Target PC and in this case 
the NSC has to be established. For this predefined SIMULINK^ blocks are available in the 
xPC Target toolbox. The communication routine used for this application is a RS-232 serial 
communication and the data was sent from the NSC to the Target PC with 2 bytes precision. 
Signals sent from the NSC to the Target PC were: 

• Mean value of the driven wheel speeds 

• Both the free rolling wheel speeds 

• Engine rotation speed 

• Eng^e torque 

When the calculations m the Target PC were done, an action was sent back from the Target 
PC to the NSC. 

[SYSTEM DESCRIPTION & MODELLING] 

One requirement of model-based control development is an accurate model of the studied 
system, which preferably is a front whed driv^ four cylinder, spark-ignition engine car. 
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Furthermore, it is of great importance to choose a model that captures the main characteristics 
of the system without becoming too large to handle in Mutations and control applications. 

[Engine] 

5 

The combustion engine converts chemical energy into mechanical energy in a four-stroke cycle 
process. The engine used in this embodiment of the invention is a turbo-charged spark-ignition 
engine also known as a gasoline engine, described in FIG, 6, Here comes a brief introduction 
to the basic principle of how the four-stroke cycle combustion works. For additional 
10 information see: 

J.B. Heywood. Internal Combustion Engine Fundamentals. McGraw-Hill International 
Editings, Singapore, 1988. 

15 The driver controls the engine torque through the accelerator that is connected to the throttle. 
The throttle controls the flow of air into the intake manifold and thus how much energy that 
can be developed in the engine. Before the air enters the cylinder it is injected with fuel to 
produce a highly inflammable mixture. The air/fiiel ratio in this mixture is carefully controlled 
to hold the so-called stoichiometric ratio. If this ratio is not reached the combustion will not be 

20 optimal and there will be an increase of emissions in the exhaust from the engine. This will lead 
to that the three-way catalyst will not work properly and it might even be damaged, depending 
on the mixture. The piston in the cylinder works in a four-stroke pattern where each stroke 
takes half an engine revolution. The strokes are called Intake, Compression, Expansion and 
Exhaust. 

25 

During intake the intake valve is open and the piston moves downward and the air/fiiel mk is 
filled into the cylinder. Then the piston moves upwards in the compression phase and the 
mixture is compressed to a higher temperature (and pressure). Just before reaching the top 
position in the cycle the mixture is ignited by a spark that initiates the combustion. The piston 
30 then moves downward (expansion) and torque is generated and transferred to the crankshaft. 
When the bottom is reached, the exhaust valve opens and while the piston moves upwards the 
combusted gas mix in the combustion chamber is pushed out into the exhaust system. Thm a 
new operating cycle starts with the intake stroke. 
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[Engine Model] 

In fhe model, the ^igine is considered as a torque generator with the ability to directly produce 
the torque requested by the driver (through the accelerator/throttle). This is not the case in 
reality where the output torque is a function of many variables such as engine speed, boost 
pressure, temperature, etc, but this it is of no concern for the realisation of this invention. 
However, the fuel-blocking method of controlling the engine torque gives rise to two 
phenomena: a rise time to reach the desired torque and an engine speed dependent delay. 

The low-pass filter step-response is due to remaining non-combusted gasoline on the wall of 
the cylinder that will ignite even though no new gasoUne is injected. Thus, 



where TdemanM is the torque requested by the controller and Te is the engine output torque. The 
time constant te of this step-response is not easy to specify but should be engine speed 
dependent to a certain extent. The conclusion of this is that even if the anti-slip system sends 
proper actions to the engine it is not certain that they can be produced in good time. 

The delay is due to the fact that the engine control system first must wait for a cylinder to 
reach the fuel injection phase before the fiiel can be blocked and then the effect of the fiiel 
blocking will take place one en^e revohttion later. This delay is referred to as the engine 
delay and the equation for the average delay becomes 



^ - ^demanded ^) 



(1) 




5-2n 5-71 
A' a, 2-0), 



(2) 



where ooe is the engine speed. The engme delay is also shown in Figure 7. The complete en^e 
model with delay and rise time becomes: 
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f.(t)=UT^(t-d^)-m) (3) 

e 

[Diivdine] 

5 The driveline is the mechanical airangement that transfers the engine torque to the driving 
wheels. It is a quite complex system that includes the dutch, a multi-shift gearbox and several 
shafts for power transfer. The car studied in this preferred embodiment of the mvention had a 
frontal drive with means that the engine torque (ideally) was evenly distributed through a so- 
called differential between the two front wheels. The main properties of the driveline consist 

10 of the gear ratios, final drive, inertias and stiffiiess and damping factors for the shafts and the 
transmis^on. 

[Driveline Model] 

15 There exist a number of possible ways to emulate the behaviour of a vehicle driveline. A three 
state model was used by Fredriksson m 

J. Fredriksson. Nonlinear Model-based Control of Automotive Powertrain. Ph.D. thesis, 
Chahners University of Technology, Goteborg, Sweden, 2002, 

20 

and it is presented in FIG. 8 in a shghtly modified form. 

It can be shown that this third-order modd captures all of the driveline major characteristics 
induding first order oscillations. Still it is necessary to assume that all other parts of the 
25 driveline (the clutch, transmission, etc.) are completdy stiff. The main difference between the 
modd used by Fredriksson and the modd used here is that the no-slip condition found in 
Fredrikssons model is replaced with a single-wheel traction model presented m the next 
section. From left to right in FIG. 8 the equations becomes 



30 
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where Je is the inertia of the engine (flywheel, camshaft, pistons, etc) and Te is the engine 
torque. For the transmission the following two expressions hold: 

CD.^io), (5) 
71^ = 2; (6) 



where i denotes the chosen gear ratio (no gear shiftmg situations are considered in this model). 
The driveshaft torque Tg, represented by a torsion spring, can be expressed as 

-k^{9s -^J (7) 

where k^s and dds are the properties of the spring and fip^ - ^ is the angular deformation of the 
driveshaft. Finally the equation for the wheel becomes, 

/>,=r,-r^^ (8) 

where /w is the longitudmal traction force from the tire model. The origin and behaviour of 
Ftrac "vvill be explained in the next section. After some manipulations the expressions become 

jA=2;-7(U<p.-^J+^*k-«J (9) 

ffl.=/ffl, (10) 



If tiie torsion of the driveshaft is introduced as a state variable as f»A = f i - f'lr the equation 
system can be rewritten as 
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(12) 



vAiGVG Te and Ftrae 316 seOT OS inputs (this is not entirely true for Ftrae as will be shown in the 
next section). Since the torque of the flexible shaft is easier to measure it would perhaps be 
better to use as state instead pf tpA but for this embodiment this is not of importance. 

[Tire-Ground Contact] 

A common simplification in mechanics is to let the peripheral speed of a rolling wheel to be 
equal to the vehicles velocity. In reality, this is true only if there is no exchange of forces 
present in the longitudinal direction between the wheel and the ground. If there are forces 
present the wheel will undergo an elastic deformation (c£ spring) and start to "creep" relative 
to the ground and so the peripheral velocity of the wheel wall differ from the vehicle speed. For 
example, to accelerate a vehicle the driving wheel speed must be higher than the vehicles 
speed. This difif^ence is called slip and is a fimdam^xtal property in vehicle dynamics. The 
SAE J670 definition of the slip A is 



where r is the wheel radii, {Bw is the speed of the driving wheel and is the vehicle speed 
(usually approximated by the speed of the non-driven wheels). When driving along a curved 
path there must be forces present in the lateral direction as well to keep the vehicle on track. In 
analogy with the case of longitudinal forces the wheel must have a lateral velocity relative to 
the groimd to transfer forces in that direction. This brings us to the definition of the slip-angle 
as illustrated in FIG. 9, vAlgtg a is the slip-angle and and Vy is the longitudinal an lateral 
velocities of the vehicle, respectively. The Mathematical expression is. 




(13) 
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a = ar( 




(14) 



As stated above both the slip X and the slip angle a are functions of the transferred force 
between the wheel and the ground. For mathematical reasons the relationship is often seen the 
other way around and a so-called slip curve can be drawn from measurements. As the graph in 
FIG. 10 shows, the friction coefficient ju is not a constant as in the case of normal friction but it 
is a non-linear fimction of the slip X (this is not entirely true and will be discussed later on). 

/I is referred to as traction coefficient to avoid confusion. The purpose of a tire 
model is to desaibe the non-linear relationship between the slip and the transferred fprce 
between the wheel and the ground. Typically the relationship can be described as 



where iw is the longitudinal force^ Ft is the present vertical load on the wheel, a is the slip 
angle (only present when comering), Vx is the vehicle speed and is the ntayimal friction 



the shape of the slip curve is dependent of the ground texture and other environmental 
characteristics, such as temperature (that influence the wheel stifBiess). However, an attempt 
to briefly e^lain why the slip curve has the form shown in HG. 10 follows. For additional 
information, see 

J.Y. Wong. Theory of Ground Vehicles, 2nd editioa John Wiley & sons. New York, 1993; 



J. C. Dixon. Tires, Suspension and Handling. 2nd edition. Society of Automotive En^eers, 
Warrendale, USA, 1996. 



(15) 



coefficient for the tire/ground combination. The word tire-model is in a way misleading smce 



and 



At low slip levels (FIG. 10, case a) there is a linear relationship between slip and traction. This 
is mamly due to the feet that the ,tire is (approximately) not sliding against the ground but 
behaves some^sdiat like a linear spring as stated above. However, the stiffiiess of this spring is 
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heavily dependent of the ground conditions, that is, the slope of the slip curve varies depending 
on \;^ether the vehicle is driving on gravel or asphalt. Alinear model for low slip values can be 
established as 

5 fi^k^,^ (16) 

vAxssQ kaip is a function of the tire characteristics (inflation pressure, tread, etc) and of 
environmental &ctors (tenq>erature, road texture, etc). 

10 As the slip X increases (Figure 4.5, case b) the relationship becomes more non-linear. This is 
explained by the &ct that the maxfanum traction (that is ;w) is reached in tail part of the 
contact patch betwem tire and ground and a sliding region will develop. 

At high slip levels (Figure 4.5, case c) the slope of the slip curve is negative. This is because 
15 the rf\^^rr\\\m traction (that is /w) is reached in the whole contact patch. When the whole tire 
is sliding heat will develop in all its length at a rate that is assumable proportional to the speed 
difierence betwem the wheel perimeter and the ground. The energy loss due to this heat 
development will bend the slip curve downwards and it assumes the characteristic humped, or 
thermal shape. However, the liump" is more or less present on different sur&ces. On low-p. 
2 0 surfaces like gravel or snow it is hardly detectable (isothermal). 

A consequence of this is that the slip curve is not only dependent of the slip X but also of the 
absolute speed difference Av. The conclusion to draw from this discussion is that the slip curve 
will have a different behaviour at different speeds since there is a sliding region developed in 
25 the contact path. 

It is also appropriate to note that most of the slip curves presented in literature are intended for 
braking applications wh^e the sUp is defined slightly different compared with the definition 
used here. This might explain the dominant "Immp" m those slip curves. While accel^ating 
30 (the case in this study) this ^lamq)" is as stated above not detected as will be shown m the 
description of HG. 20. 
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[Tire-Ground Contact Model] 

The basic behaviour of a wheel driven vehicle travelling on a straight line is quite well captured 
by a single wheel or quarter car model (FIG. 1 1) as described in 

5 

T. Glad, L. Ljung. Reglerteori. Studentlitteratur, Lund, Sweden, 1997. 
The equations of the model are (cf equation 8): 



radii, m is the vehicle mass, is the vehicle speed, iw is the force in the contact path, Fr^ is 
15 the rolling resistance, Q is the air resistance (e.g. drag) coefficient, A is the vehicle frontal area 
and /> is the air density. 

The longitudinal traction force Ftne could as stated earher be seen as a non-linear ftmction of 
the slip with a number of parameters^ such as temperature and ground sur&ce conditions. One 
20 approach to the use of the slip curve in a simulation sense is to say that the shape of the curve 
is depending of the environment (temperature, ground surface conditions, etc) and then use the 
slip X to calculate the longitudinal force. 

The non-linear relationship between the longitudinal force and the slip is then preferably 
25 calculated using a very simple static tire moddi as 



10 



(16) 



(17) 



where is the wheel inertia, q)^ is the wheel speed, is the applied torque, Tw is the wheel 




if X^X^ 



(18) 



(19) 



30 
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The parameters of the equation could then be changed to simulate drivmg on diflferent surfeces 
(generally a sur&ce with lowA/^ has a low .^). 

[Complete Model] 

The resulting vehicle model consists, even though kept simple, of a five-state system of 
diflferential algebraic equations. The model was implemented in SIMULINK® and served as a 
primary environment for the development of the control system of the preferred embodiment of 
the mvention. Tdestnd is the input to the system. 






(21) 




(22) 



(23) 



[DESIGN OF CONTROLLERS] 
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The quality of the input signals to the controller is often a limiter when it comes to how refined 
the controller is possible to make. For instance, many controllers make use of the derivate of 
the controlled entity to be able to fester compensate for changes. Perhaps even fester 
compensation could be reached if the second derivate also was used but this is seldom possible 
due to the feet that differentiated agnals have a tenden«y to become very noisy. 

[Model simplification] 

To derive a manageable model-based controller it is appropriate to simplify the equations for 
the considered system (here, equations 20-23). This is a hard trade off between keeping 
information and still getting a nfianageable system. -Two different embodiments for the 
controller are here proposed. They both use the same simplifications (except for the slow 
variation of vehicle speed which is only used by the normalised controller below): 

• The driveshaft is stiff 

• Engine dynamics can be sufBciently captured by a LP-filter and a delay 

• Linear tire model up to /w 

• Slow variation of vehicle speed 

[Stiff drive shaft] 

Test performed on the V40 shows that the driveshaft is quite stiff and the wheel oscillations 
can be ignored. So for the controller design purpose it is not necessary to include the influence 
firom the **weal^' driveshaft. 

This makes it possible to transform the inertia from the engine to the wheel over the 
transmission leading to a new inertia J that describes both the two wheels and the engine 
inertia as 



7 = Z/,+/V,. 



(24) 
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\^th a stiff drive shaft the differential equation for the driveline (equation 12) is derived as 



The stiff driveshaft model is schematically viewed in FIG. 12. 
[Enguie dynamics] 

The engine dynamics due to wall-wetting and other thermodynamic phenomenon are assimied 
to be of LP-character but are very difficult to estimate. Wall wetting arises since some gasoline 
always sticks on the inlet walls to the blocked cylinder. This dynamic behaviour is actually 
heavily depending on engine speed, wall temperature and which control action being applied. 
Therefore an average estimate of the time constant 4 is preferably used. The equation is the 
same as used in the simulation environment 

[Slow variations of vehicle velocity] 

The slip function (defined in equation 13) is non-linear, therefore a direct application of a linear 
controller is not possible without some kind of simplification. In the later on described absolute 
controller this simpUfication is not necessary. The idea is to use the feet that changes in vehicle 
speed is much slower than the changes in wheel speed for the driven wheels. 

The validity of this statement is also relevant for this embodiment, although the engine inertia 
Je makes the total inertia larger than in a braking manoeuvre. 

[Linear tire model] 

In comparison to theory, the presumption of a linear tire model (FIG. 12) is a major 
simplification. On the other hand, measurments performed by the inventors imply that the 
sdnoplification may be quite relevant (see FIG. 20). Anyway, in reality the sUp curve is not fixed 
but changes all the time (rolling gravel on asphalt, tram tracks, etc.) so even a non-Unear tire 
model would not give a very good picture of the reaUty. A more sophisticated strategy for 




(25) 
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estimating the slip curve is disclosed in an alternative embodiment which follows later on. 
However, even if the linear slip curve isn't entirely correct, it provides usefiil qualitative 
information about the groimd conditions. 

5 In a possible embodiment of the invention two controllers are used, one controller on the top 
of the tire model \diere 



^ = 0 (26) 



10 and another controller on the slope where 



•^ = *^ . (27) 



15 



20 



and 

Fir^^J^supF^ (28) 

FIG. 13 shows the relationship between the real slip curve and the approximation used in the 
controller. 

[Normalised slip transfer function] 



In the transfer fimction that is to be derived, the common way to define the slip A is used, that 
is the normalised slip definition (equation 13). The assumption that the change in vehicle speed 
25 is much slower than the change in the slip gives that the vehicle acceleration can be seen as 
zero. The derivative of the slip then becomes 



A = ^2£^k. (29) 



4 
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Combining equations 3, 25, 28, and 29 and transfonning them to the Laplace plane, leads to a 
transfer function from demanded engine torque to slip as 



_ 5x_ 

X ie 
G^(s)=^^;r^ = -r-^ : 



(30) 



Unfortunately this equation has four time varymg parameters, the vehicle velocity V;^ slip slope 
ksup, engine speed Oe and gear ratio i. For simplicity, the engine speed can be seen as a function 
of the vdiide speed and gear. This is of course not entirely true, especially in a high slip 
situation but this assumption gives some extra margjns in the controller when it overestimates 
10 the engine dday. For the other three parameters we need to scale the system. This gives us a 
gain scheduling controller in the three dunensions vehicle speed, slip slope and gear. The slip 
slope is preferably estimated with a Kahnan filter, as will be described later. 

The transfer flinction Gj^^s) was compared with the behaviour of the real vehicle and the 
15 result is diagrammed in FIG. 14. The transfer fimction was fed with the measured engine 
torque, vehicle velocity and the estimate of the prevailing ksiip. 

[Absolute slip transfer function] 

2 0 Another way to avoid non-linearity m the controller is to use the absolute slip definition, that is 
the speed difference iiv between the diivm and the non-^ven wheels, 

Av = a?^^-v^. (31) 
25 The relations betwem the traction ju and the speed difference is 

;i = *A.Av. (32) 
wh^e kAv can be derived from the estimat 

30 
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(33) 



Equation 31 is differentiated in analogy with the procedure for the normalised slip transfer 
function although it is not necessary to assume that the vehicle speed is constant. Thus, 



(34) 



10 



Togeth^ with simplified fimctions for the driven wheel speed (25) and the vehicle velocity (cf. 
equation 20), 



m 



(35) 
(36) 



a transfer function for the absolute slip can be derived by combining equations 3, 23, 32, 33, 
15 34, 3 S and 36 and transforming them to the Laplace plane. The transfer fimction becomes 



Av 



V. 



Av 



(37) 



with the time varying parameters i and hi,. The transfer function Ga,^s) was conq)ared with 
20 the bdiaviour of the real vehicle and the result is diagrammed in FIG. IS. The transfer fimction 
was fed witli the measured engme torque and the prevailing k^ A direct conqpaiison between 
the two transfer fimctions and thdr ability to reproduce reality (FIG. 14 and IS) implies that 
tiie transfer functions for the absolute slip seems more accurate. This is probably due to the 
feet that the normalised transfer function looses information when the vehide speed derivate is 
25 omitted. 



[Controller design - IMC] 
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The method preferably chosen for designing the controDer in hereinafter referred to as Internal 
Model Control (IMC). When choosing a design method it is important to consider the problem 
specific criteria that have to be met. For this application three main criteria could be 
established. 

First, the controller should preferably be quick but still robust and preferably be able to follow 
and hold a reference vahie. Generally this can be achieved by using a controller that includes 
parts proportional to the error (P), the integral of the error 00 and the derivative of the error 
(D) or equivalent. 

The second criterion arises since the system (the vehicle) will change its behaviour depending 
on the working point (that is gear ratio, groimd conditions, etc.). This implies that the 
controller should preferably have to be gain scheduled to obtain a good performance. In order 
for this to work properly, the derived controller should preferably be "smooth" (or even linear) 
in its components with regard to its scaling parameters. This is to avoid sudden changes in the 
controller paramet^s when changing operating point with resulting "jump/' behaviour of the 
control signal. 

In addition to this a third criterion was that the chosen method should preferably be easy to use 
and result in an easily tuned controller. 

The general concept in IMC is to design a controller that gives the closed loop system (FIG. 
16) a desired transfer fimction. If the transfer fimction, describing the controlled system G(s) is 
completely correct this method automatically results in a stable controller with regard to the 

sensitivity function and the complementary sensitivity fimction (Sensitivity fimction S(s) = — , 

V * 

Complementary sensitivity fimction Q(j) =— ). The controller has, in this embodiment, 

w 

preferably only one design paramet^ fi that sets the bandwidth for the system. 



wo 2004/067307 



-29- 



PCT/EP2003/000807 



To illustrate the influence of P step responses for a closed loop system is diagramted in FIG 17. 
The graphs show that there exists a that results in the most &vourable response. For more 
information about IMC, see 

T. Glad, L. Ljung. Reglerteori. Studeotlitt^tur, Lund, Sweden, 1997. 



The following controllers, developed with IMC, can be seen as PID-controIlers convolved with 
first order low pass filters. Except for the design parameter all parameters in the controllers 
have physical meaning describing the slip dynamics. It is also clearly visible how the different 
model parameters and a affects the P-, I- and D-parts of the controller. For the normalised slip 
the controller becomes 



Tifir^ 



5+ 



(38) 



which needs to be scheduled for v, ka^ and t For the absohite slip the controller becomes 



Jt. 



5 + 



F,8 



^rl 1 
— +— 
J m 



JJl L 



AtrJ 



s 



(39) 



which needs to be scheduled for Ar^y and /. 



The effects of the engine delay (and communication related delays) are not treated in any 
special way, but the controller is slowed down to keep stability. There are other ways to handle 
time ddays, the most common is to use a prediction model as in the Otto-Smith controller 
described in the reference dted above. 



[ControUer Logic] 
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The main logical problem for an anti-slip system is how and when the controller should be 
turned on and off. The first subquestion to this problem is what the controller should calculate, 
the necessary torque the engme should give to sustain a certain slip or the torque reduction 
necessary in the current situation. Another related question is how the integrator should act 
^en the controller is outside its working area. 

[Ackerman Steering] 

When driving along a curved path, the front wheels of the vehicle will travel a longer distance 
than the rear wheels. If the car is front wheel driven this will increase the slip. To compensate 
for this "false" slip component it is advantageous to recalculate the shp with regard to the turn 
radius. If the slip angles are neglected (in reality they are preferably kept small) and no driving 
torque is applied (Vfiont = eWw) a bicycle model will turn like in FIG. 18 and the velocities in 
the slip equation could be normalised with respect to the travelled distance for each tire. 

This compensation yields approximately even for small slip and slip angles. The compensated 
slip equation then becomes 



r CO 



r CD 

W W 



(40) 



where L is the shaft distance and R is the turning radius. 

This compensation is dependent of the unknown turning radius R so it has to be estimated. One 
possible way of doing this in a vehicle not equipped with accelerometers and yaw rate sensors 
is to derive R from the difference in rear wheel vdocities as 

where ^ is the trade width of the vehicle. The conq)ensated slip equation can be established 
like 
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(42) 



If absolute slip is used the compensation can be derived in analogy with the discussion above 
5 and the compensated Av becomes: 



The line of thought used to derive this slip compensation is quit simplified compared to reality 
10 where both slip angles, and speed differences between the firont wheels are present. The 
cornering behaviour of a real car is also dependent of where its centre of mass is situated as 
also described in 

J.Y. Wong. Theory of Ground Vehicles, 2^ edition. John Wiley & sons. New York, 1993. 

15 

[The ksiip-estimator] 

The only unknown parameter in the controller is the slip slope k^^ that has to be estimated. 
This is not a simple task to perfonn in real time applications but there are some methods 
20 developed for this task, an example is RFI (Road Friction Indicator). The RFI is preferably 
based on a time varying Kalman filter as described in more detail later on. In the present 
embodiment, a simplified scalar version of the Kalman filter is used. 

The filter is advantageously based on the following equations where eft) and vft) are 
25 measur^ent and process noise: 



Av = 




(43) 



(44) 



A=A^(0-V(0+KO 



(45) 
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Q(t) and R(t) is introduced as the variance of v(t) and e(t): 



Q = 4emt)] (46) 
/? = E[v(Ov(0] (47) 

>^th v(t) and e(t) considered as white gaussian noises the Kabnan filter gives the optimal state 
estimate k(t) . The filter algorithm advantageously looks as follows: 

S(t) = P(t-l)+Q (48) 

m=s(t)M(tMm)M(t)+Ry (49) 

^^,(0"' =k^(t-ir +K(t^(t)-Mit)k^itr) (50) 
Pit) = Sit)-K(t)Mit)Sit) (51) 



The design parameters for the filter is the ratio of the noise cx)variances Q and R (the initiations 
values oiP(0) and ksiip(0) can be stored when the car is not running). The covariance ratio is 
set so that k^i^ is quite rapidly estimated but still not too noisy. For sudden changes in k^ip (e.g. 
ice patch) the estimator will not be able to capture the rq>id change and thus the controller will 
not be properly initialised (the k,i^ used in the controQer and the slip reference calculator will 
be wrong). 

As for the controller, the Aj/^-estimator should preferably be surrounded by logic to switch the 
filter on and off. 

[REFERENCE VALUE CALCULATION] 

In classic theory the slope kaip of the slip curve is said to depend only of the tire properties 
(essentially inflation pressure and tire type). This has bem taken for a fact for many but the 
mventors have found out that this is not at all the case. On the contrary, the slope is heavily 
dependent of the ground conditions. When driving on gravel or ice the estimated kup is found 
to be substantially lower than what estimated while driving on asphalt. This is deariy 
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contradictoiy with the commonly accepted theories on the subject. With this new theory, an 
on-line estimate of a ground condition dependent kji^ could, together v/ith empirical 
kno\'^edge of the corresponding fi^, be used to estimate an appropriate slip reference for the 
current driving case. This approach is preferably used in this embodimoot because of its 
simplicity even though an alternative embodiment is discussed later oa 

[Implemented Strategy] 



The slip reference is advantageously calculated according to two assumptions: 

• The incline of the slip curve can be used to set a reference value (empirical knowledge 
of the corresponding ;w) 

• Lateral stability is prioritised while cornering and at high speeds 

The calculations are preferably performed as following: first a nominal reference value is 
chosen dependmg on prevailing ground conditions (e.g. kup) as 



\jaw if ^tHp ~ ^sttpjaw 



(52) 



with some kind of linear interpolation betwem the entries. The required lateral force to ke^ 
the vehide on the path is calculated and the reference vahie could be conq)aisated as 



(53) 



According to equation (41) R can be derived from the wheel speed dgnals and the 
compensated reference value becomes 



(54) 
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vAxero kj and A2 are the design parameters. Even if the safety ellipse is a good model for 
deriving the optimal reference value it is dependent of many parameters that are hard to 
estimate. The equation suggested here is more intuitive, the reference value will decrease with 
larger centripetal acceleration but it will never go down to zero (and thus kill the engine). 
However, there is no possibility of detecting what the driver intends to do in the sense that the 
steer angle is not measured and thus it is impossible to detect a turn imtil a certain speed 
difference is reached on the rear wheels. This implies that if the driver intends to turn he has to 
go easy on the accelerator to not provoke large slip angles on the front wheels. If there are 
large shp angles on the front wheels the available lateral force will not be sufficient to turn the 
car and no reduction of the reference value is performed with the result that an uncontrolled 
slide is provoked. 

Finally, if desired, the speed influence could be considered and the reference value 
compensated according to 



AvM (55) 



where ks is used to control the speed influence. For low speeds anti-slip system should 
preferably work as a traction system, that is, the slip reference should preferably be set high 
while side stability is prioritised at higher speeds. Altogether this formula becomes very easy to 
understand and the influence of the different compensations can easily be tuned with k], kg and 
ks to achieve the desired behaviour of the vehide without any knowledge of the theoretical slip 
characteristics. There is also the possibility to let k and Aj depend on kup if there are 
dmiands for, let's say, even better lateral stability \>^en driving on ice. 

The calculated normalised slip reference could then easily be transformed into a absolute slip 
reference by multiplication with the vehicle speed. 

[Adaptive slip reference calculation] 
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The embodiment suggested above does work sufficient for one set of tires or tires with similar 
characteristics, but the reference vahie may become incorrect if the tire type is radically 
changed (e.g. switch between summer and winter tires). 

Another, simpler, ^proach to the problem with dififerent tires could be to derive some 
standard settings for summer and winter tires and then let the driver decide his preference by 
pressing a buttoa 

[IMPLEMENTATION AND PERFORMANCE] 

This section will deal with the aspects of implementation of the control system in the Volvo 
V40 T4. This vehicle is originally equipped with an anti-slip system (NSC) which is based on 
look-up tables. 

The developed control system could be directly implemented in any frontal driven vehicle with 
fuel blocking possibility, although it has to be adapted m the sense that the car spedfic 
parameters such as mass, wheelbase, inertias, etc., have to be changed. These parameters are 
usuaUy found in the car configuration files sent over the CAN bus. 

[General Implementation Aspects] 

[Delays and Sample Time] 

The data processing and data transfer will add an extra delay time to the control loop. In total 
the delay will consist of two parts, a constant part due to the vehicles intamal signal routing 
configuration and sensors and a varying part that will depend on the engme speed (de^^te as 
described h^einafter within the scope of the engine model). Thus, * 

dtai=d^-^dengine(Pe) (56) 

The lengtii of the constant part of this delay is difficult to estimate when it is dependent of how 
the different subsystems work. It could be estimated by triggering the ABS-system and then 
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measure the lime elapsed before the EECU delivers an output to the engine. However, when 
developing the system the delay will be significantly larger compared with the delay of a 
product unplementation due to the use of xPC Target. 

5 In addition to the extra delay, xPC Target cannot work faster than 50 Hz. This is not due to 
xPC Target itself but to the feet that the NSC cannot transfer data fester than this rate. The 
existing anti-slip controller (NSC) works at 200Hz. This implies that a direct comparison 
between the two systems is not completely fair. 

1 0 [Parameter Identification] 

The model was validated against real-world data in a simple experiment. Since the 
measurement equipment was not able to preset negative torques firom the engine, i.e. all cases 
of engine brakmg, only acceleration cases could be validated with any certainty. Zeros replaced 

15 negative torques. Since influmce form ground surface conditions were undesired the fi'ont 
wheels of the V40 was elevated in mid-air during testing. Full throttie accelerations were 
performed and the results presented in FIG. 19 as solid lines. The model was then fed with the 
measured torque output firom the real engine and the model output was diagrammed against 
the measured outputs. The model parameters were adjusted fi^om their (estimated), physical 

2 0 values to better reproduce the behaviour of the vehicle. 

The most important model parameters to estimate were J« and as they are crucial for the 
controller performance. As dearly viable m FIG. 19 the model diverted severely fi-om the real 
vehicle when the throttie was closed (negative engine torque). 

25 

[The Measured Slip Curve] 

When confironted with reality one discovers that the tire model ^ves a idealised version of the 
groimd contact characteristics. The measured slip may contam a lot of high fi-equency 
30 components due to the driveline dynamics, oscillations of the tire body and variations in the 
ground surfece (stones, bumps, holes, etc). The longitudinal traction force is not directiy 
measurable but has to be derived firom either the vehicle acceleration or the engine torque. 
Both methods have their drawbacks. The acceleration is preferably measured as the derivative 
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of the rear wheel speed and thus quite noisy. Furthermore, the method according to the present 
embodiment does not c^ture influence from slopes, air resistance, etc. The engine torque from 
the EECU is m the preferred embodiment taken from a steady state and cannot capture 
any transient behaviour of the engine. In addition to this the engine torque is only valid for slip 
curve estimation at low slip levels if it is not correctly compensated for the wheel inertia and 
other energy losses. In FIG. 20 the vehicle acceleration was used as input. Both the slip and 
the acceleration had to be LP-filtered to be able to draw the slip curve. 

However, the resulting slip curve shows that the assumption of a linear tire model is quite 
accurate up tiUl Urnax and that no significant "hump" could be detected. Thus the strategy with 
two controllers, one in the linear region and one in the unstable, as advantageously suggested 
by the present embodiment seem promi^g. 

[kdip-estimator] 

The filter is preferably tuned to have a feirly rapid transient behaviour without giving a too 
noisy estimate of the slip slope. As shown in FIG. 21 it takes about four seconds to detect a 
change in ground condition (here asphalt to gravel). However, additional work has to be done 
on the estimator if it is to be more accurate. The flats in the curve indicate where the filter logic 
has decided to turn the estimator off. 

[PEKFORMANCE OF CONTROIXER] 

Of the two described embodiments for the controllers, the controller based on absolute slip 
seemed to be most suitable for vehicle implementation mainly due to better results in the 
sunulation environment. Di addition to this, by using absolute slip 2lv a couple of other practical 
problems could be handled in a smoother way. For instance, the nonnalised slip becomes 
undefined when the vehicle speed is zero and very noisy at low speeds. 

The controller was bench marketed against a DSA, the existing anti-slip system implemented in 
the Volvo SA^40-series. 



[Full Throttle Takeoff] 
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From an initial engine speed of 3000 RPM, a hard take-off was performed at full throttle in 
first gear. The results are presented in HG. 22 and FIG. 23 for DSA and FIG. 24andFIG.25 
for NSC 2. 

5 

[Acceleration] 

The car was slowly brought to an initial speed of 30 km/h in second gear. The data logging 
was started and fiill throttle acceleration was performed up to 60 km/h. The results are 
10 presented in FIG. 26 and HG. 27 for DSA and FIG. 28 and FIG. 29 for NSC 2. 

[Cornering] 

The car was slowly brought to an initial speed of 20 km/h in first gear. The accelerator was 
15 fiilly released and the steering wheel was turned to its maximum angle. Then the data logging 
was started and the accd^ator was pressed to fidl throttle while the car completed two turns 
in a circular path. After the two complete turns, the car was straigjitened up and accelerated at 
full throttle to an engine speed of --5000 RPM. The results are presented in FIG. 30 and FIG. 
3 1 for DSA and FIG. 32 and FIG, 33 for NSC 2. 

20 

[Acceleration Lane] 

An attempt was made to imitate a more realistic driving scenario when the vehicle is 
accelerating through a long curve, sudi as an acceleration lane for entering a highway. The 
25 vehicle was driven at fiiU throttle in third gear. This test was difficult to repeat in exactly the 
same manner for the two controllers so the two test runs are not completely comparable. 
The results are presented in FIG. 34 for DSA and FIG. 3 5 for NSC 2. 

[Summary of test results] 

30 

When evaluating the performance of anti-slip systems there a advantageously separate issues to 
consider, comprising how the actual slip behaves and how the driver experi^ces the vehicles 
behaviour. 
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From the measured test results it seems that the controller according to the embodunent 
maintains a smooth slip vahie with less oscillating behaviour than DSA. In addition to this, the 
control signal also shows a smoother behaviour, which is preferable if the controller should 
have to use another torque reduction method such as throttle- and/or ignition angle control. 

However, the real test for anti-slip system is how the driver experiences them. When driving, 
DSA is experienced as slightly smother than the presented controller. This is probably due to 
the fact that DSA uses a higher slip reference value than the presented controller. At high slip 
levels, the difference in vehicle acceleration is small for oscillating slip while it become larger 
when the slip is kept at the slip slope. However, a high reference value tends to prioritise 
traction ov^ side stability and also increase wear on the tires so here a trade-off has to be 
made. 

[Second System Embodiment] 

Note, that in the embodiment described hereinafter, the slip /I is denoted by .y and the transfer 
function G(s) by G(p). 

[Outline of the controDer structure] 

For the sake of sunplicily, we will here repeat the basic dynamic equations assuming a one- 
wheeled vehicle. The vehide specific details for cars and motorcycles will be presented 
herdbnafier 

Assume that this vehicle has mass /x, normal force and the wheel has moment of inertia J 
and radms r. The basic dynandc equations of motion are: 
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Here fi^s) is the nonnalised traction force F/F^, The slip definition is as given in the above 
embodiment (see equation 13). This static rdation works fine fi^r our purposes. The 
knowledge of the fi relation is critical as noted above. The prindpal form of this relation is 
shown m FIG. 10 and the best possible fiiction is attained for the slip yielding tiie peak value of 
this curve. The basic control strategy is to reduce the engine torque and thus /ik so as to never 
pass the peak value. 



This gives the dynaniics for the slip as 



s = — OJ -Vj. 



That is, each tire is modelled as a first order non-linear ordinary differential equation. 
Standard control techniques apply on linear differential equations, and the approach that will be 
taken is to use a robust controller on the linearised dynamics. This involves computing a linear 
relation /i^s) + fj^^s^ (s - ^o) around the current working point Sq. 



A diivelme model fi-om the engine torque reduction signal to the tire force Ttaction similar to the 
above embodiment can be taken as 



where G(p) is a linear transfer fimction as described in detail in the above embodiment, here 
expressed m tiie differentiating operator/?. We will assume that a torque reduction algorithm is 
given. 



Tcagi«. = min(K,l)Tcommandcd, H > 0, 

where w is tiie reduction fector (could be discrete). The minimum operator is included so that 
the controller can be active all the time and normally strives to increase the engine torque (u > 
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1) SO that the slip would approach its reference value. Common principles for torque reduction 
include controlling the throttle or ignition angle or to use fiid blocking. 

A part of the algorithm is a &st model-based £riction estimator. Other intended features is a 
cascade control similarly to the above embodiment wh^e the inner loop controls each 
individual wheel slip to a reference value, and the outer loop determines a suitable reference 
value yielding a constant safety margin during cornering. 

An alternative structure of the controller, as illustrated in FIG. 36 is as follows: 

• An optional friction detector/estimator is run all the time. It provides an estimate of the 
maximal utilizable friction force /w and the initial slope k of the slip curve. This 
algorithm will be referred to as RFI (Road Friction Sidicator). 

• Whenev^ the traction (or braking) force gives a slip larger than a threshold value, the 
slip curve estimator is enabled, preferably using /w and k from RFI for initialisation, A 
good initialisation prevents divergence problems and speeds up convergence. 

• The friction ellipse describes how the total friction force can be split up between lateral 
and longitudinal force, see FIG. 37. Here it is a circle with radius 1 indicating perfect 
friction, but for slippery roads and also considering a certain safety margin, the radius 
must be decreased accordingly. For a given lateral acceleration fiy^ the maximum 
longitudinal force given the friction and safety margin, can be computed. For 
MASC, it is used to compute how much of the utilizable force that is left for traction, 
taking the driver commanded ste^g into account. The lateral acceleration is taken 
either from an accelerometer or computed from the steering wheel angle or the wheel 
speed signal, or finally using a combination of several principles in a Kabnan filter. The 
output is a suitable longitudinal reference value of longitudinal normaUsed traction 
fi>rce, 
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• An on-line estimated fiiction curve fi^^s) is estimated. Here the peak value is most 
important. From this, the reference longitudinal force is mapped to a reference 
slip j^. 

• A closed loop controller change computes a suitable engine torque to make each wheel 
slip equal to the reference value. The control signal u can be thought of as the engine 
torque reduction factor, since the commanded torque will never be increased of course. 
The reason for computing u> 0 and then Umiting it to min(i/,l) is to be able to use 
standard linear control design and to avoid non-linear switches inside the controller. 

[Non linear friction model] 

Crudal for the performance is the knowledge of the jurs relation, defined by the slip curve 
fij[sy We will start with defining useful properties of the model structure that is used in the 
estimation, and then examine some adaptive filtars and their performance on simulated and real 
data. We assume that the model structure can be expressed in the form Ms;0), where 0 
denotes the firee parameters. 

[Model requiremmts] 

Before listing wanted properties on the model structure, some reflections on the curves in PIG. 
10 are appropriate. The most characteristic properties for such curves are: 

• The initial slope k = d^J(sy{b, for.s = 0. 

In our case. A: can be assumed a priori from the RFI algorithm, while in the clasaical 
interpretation, it is a tire property not changing much over time. 

• The maximum utilizable fiiction force /w = max« ^(s). From the RFI algorithm we 
have an initial value for estimation as either high = 0.7, say) or low = 0,2, 
say). 
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• The slip for which the maximum utilizable fiiction force is obtained, Sjoa arg max. 
fijJis). A common value of ^nua in practice is 0. IS. 

• The stationary utilizable fiiction force 

These are perhaps the most important characteristics of the pjj(s) curve. Note that feasible 
values of these parameters are given by the open box 

10 

0< fc, 

0 < flrnax <h 
0 < Sfnax < OO, 
0 < Moo < 1- 

For other vahies, the corresponding curve does not correspond to physical prior knowledge of 
fiictioa The vector /WAiafc Ao) will be refisrred to as the physiccd parameters. We will 
15 look for model structures which have four parameters in 0, so we can interpret the result in 
terms of the physical parameters. There are many advantages of this, which will be pointed out. 

The most important properties on a model structure to be a good candidate for on-line 
estimation are the following ones: 

20 

PI Flexibility. A large class of empirical relations should be able to be approximated well 
by adjusting the parameter values. 

P2 The parameters 0 should be sinq)le to relate to the following physical parameters. 
25 There are two motives for a simple relation between 0 and (k, /^o): 

P2a Initialisation of the adaptive algorithm is greatly simplified if there is a simple 
relation 0 ^JQc /w, ^Wx): 
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P2b Feasibility projection. The function must of course correspond to a 

feasible slip curve. If Aere is an explicit transformation (A; /ibuuo s^gu Ao) = 
^6), then it is simple to constrain the parameters so that the physical 
parameters belong to the open box defined above. Furthermore, it simplifies 
if the feasibility region for ^has a simple form, preferably without non-linear 
couplings between the parameters. 

P3 Differentiability with respect to s. The adaptive slip controU^ is a function of 
dnJis)lds^ which thus must be computed in the controller design. 

P4 Enablmg effident identification. The best model structure in this respect is one 
which is linear in the parameters. The next best model structure is one which is 
diff^entiable in the parameters enabling gradient based adaptive algorithms. For 
gradient based algorithms, the Hes^an should be as diagonal as possible, which 
practically means that the parameters should be as decoupled as possible. 

[CARMA model] 

We here examine a structure which is linear in the parameters. The step response of a second 
order Continuous time Auto-Regres^ve Moving Average (CARMA) model, as often used in 
control theory, has a sunilar shape as the slip curve. Think of slip as time' here. In the 'time' 
domain, the step response is 

= (A; + k)sx{l - e^^/'^) - jksaCl - e^''^^). 

A non-minimum phase zero causes an overshoot in the step response, which models the peak 
of the slip curve. One advantage of using model structures borrowed firom control theory, is 
that standard system identification algorithms may be used. For instance, the pern function of 
the system identification toolbox might be modified to "non-causal" sampling. This particular 
parametrisation gives (P3) 
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from which we compute the physical parameters 

A*'(0) = * 
tix{oo) = {k + k)si - k82 



This means that P2b is complied. Unfortunatdy, the feasibility region is hard to determine aid 
there does not seem to be an explidt inverse, so P2a is not satisfied. The parameter gradioit 
P4 has a quite simple form, enabling computationally sinq>le adq>ttve algorithms. 



(Tiie peak parametrised model] 

The second proposal is directly parametrised in the peak parameters (fj„^ s„^), 

Ma)=/w.ri-i4^-iry 

The expression inside the absolute sign is positive when s > and negative otherwise. 



Assume first, for simplicity, that n is constrained to be even. The derivative is then given by 



The physical parameters can be conq)uted as 
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Note that the straight line between the origin and the peak has slope fimJs^, The factor 
n(l^r)/r is thus the excesaive slope, and by feasibility this must be larger than one so we must 
have n > r/{l+ry From the expression for we interpret as the peak excess, compared to 
the stationary value. We must have r ^ 1, so the feasibility region is a unit cube for 
r). To sunomarize, we should design the adaptive filter to guarantee that 

0 < r < 1, 
0 < AWc < 1, 
0<s™„<l, 

Q<n< — 

l + r 



One drawback with this model is the non-linear coupling between (/i,r) and which are 
influencing different parts of the slip curve. That is, property P4 might not be perfectly satisfied 
here. 



[Adaptive slip curve estimation] 
[Adaptive algorithm] 

The peiiiaps simplest possible adaptive algorithm for the peak model is a stochastic gradient 
algorithm: 



e=et 



wo 2004/067307 



-47- 



PCT/EP2003/000807 



Here a is a step size to be chosen. The algorithm should only be run when the sUp is larger 
than a certain threshold value, say 5 > 0.05, to avoid excitation problem. For small slip values, 
only the slope is identifiable. At each iteration, the update is accepted only if each updated 
parameter belongs to the feasibility region. Instead of computing the analytical gradient, it 
might sunplify to use a numerical gradient. 

[Simulation results] 

The performance of the numerical gradient algorithm is here investigated on simulated data. 
The system is assumed to be in closed loop, so the slip is cyclically swept over the peak value, 
here according to a sine wave after an initial ramp, as shown in the first diagram in FIG. 39. 
The second diagram shows how an initial guess converges to the true curve, and the last 
diagram shows how the parameters converge. 

FIG. 40 shows how the algorithm handles different initial conditions. Most important is to 
compare the second last diagram to the final estimated slip curves in the last diagram. 

The left diagram of FIG. 41 shows the trajectory for s and ju on left and right driven wheel, 
respectively, during a drive with fiiUy developed spin. First subdiagram for static data, and 
second subdiagram for LP-filtered fi. Right diagrams show time diagrams of s and /i and 
estimated slip curve. Three different test drives are compared in the different rows. 

[Other Embodiments] 

Here some other embodiments of aspects of the mvention are presented. 
[Controller using internal states] 

If absolute slip Av is used and air resistance is neglected it is possible to set up a linear model 
for the vehicle with the states xj = Te, Xf^cne, x3 - (o^x4 ^A(p and x5 = v,. A Hy- or Ho- 
controller could be based on this system. However, such a controller would demand a quite 
extensive estimator (e.g. Kahnan filter) to keep track over the non-measurable states, Te, and 
A(p. The vehicle load could then be estimated. 
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x = Ax+Bu+Mv with K = r„g and v==/5to„rf (57) 
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[RPM Control] 



One problem in slip control is the influence from the dynamic properties of the driveline. The 
problem consists of the fact that the controlled entity (the slip) is situated away" from the 

10 actuator (the eng^le) and oscillations from the driveline will have a negative effect on the 
stability of the system if not correctly modelled. In addition to this the slip-signal will be quite 
noisy while the RPM signal is quite "clean" in comparison. If the engine speed (ELPM) was 
controlled instead of the wheel speed those problems could be avoided. However, other 
problems arise such as the reference value calculation that will become highly gear ratio 

15 dependent and how to handle driving cases with active use of the clutch. 



[MISCELLANEOUS] 



According to an exemplified embodiment of the invention, the modd based anti-slip system 
20 may be any machine cq)able of executing a sequence of instructions that specify actions to be 
taken by that machine for cau^g the machine to perform any one of the methodologies 
discussed above. The machine may be an application specified integrated drcuit (ASIC) 
including a processor and a memory unit. The instructions may reside, completely or at least 
partially, within the memory and/or within the processor. 
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In particular, the model based anti-slip system may be implemented in the form of a computer 
system within which a sequence of instructions may be executed. The computer system may 
then further include a video display unit, an alpha-numeric input device (e.g. a keyboard), a 
cursor control device (e.g. a mouse), a disk drive unit. The disk drive unit includes a machine- 
5 readable medmm on which is stored the sequence of instnxctions (i.e., a computer program or 
software) embodying any one, or all, of the methodologies described above. 

The computer program product may be a machine-readable medium which is capable of storing 
or encoding the sequence of instructions for execution by the machine and that cause the 
10 machine to perform any one of the methodologies of the present invention. The machine- 
readable medium shall accordmgly be taken to mdude, but not be limited to, solid-state 
memories^ optical and magnetic disks, and carrier wave signals. 

In the above description, for simplification, it is mostly referred to the method case. The 
15 system case can be easily derived fi-om the latter one by replacing at the appropriate places in 
the description the expression *step' for the method case by the e?q)ression 'component* for the 
system case. 

Furthermore, all other publications and existing systems mentioned in this specification are 
2 0 herein incorporated by reference. 

Although certain methods^ systmis and products constructed in accordance with the teachings 
of the invention have been described herdn, the scope of coverage of this patent is not limited 
thereto. Particulariy, the Volvo V40 T4 environment serves solely as an exemplifymg 
25 embodiment of the present invention. It is explicitly pointed out that the invention is in no 
respect limited to this particular vehicle environment. On the contrary, this patent cov^ all 
embodiments of the teachings of the invention fiiirly falling within the scope of the appended 
claims. 
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CLAIMS 



1. A method for coxrtroUing the wheel slip (X) of at least one wheel on a vehicle, 
comprising the steps of: 

- obtaining at least one dynandcal state parameter of the vehicle, 

- determining (G(s)) a slip reference value indicative of the wheel slip at a determined 
driveline torque on the basis of a dynamic model of the relevant vehicle 
components, and 

- controlling (R(s)) the driveline torque (L) of the vehicle by determining a torque 
manipulation signal based on the dynamical state parameters and the sUp reference 
value. 

2. The method of claim 1, \^^erdn the model-based calculation is based on a vehicle 
model, comprisang at least one of 

- an engine model, 

- a dynamic model of the driveline for estimating forces in the tires' contact P^th^and^ 

- a tire-road contact model. 

3. The method of claim 2, wherein that the tire-road contact model comprises a slip- 
friction force model resulting in a slip curve indicative of the respi^venroad 



5. The method of claim 3, wherein the slip-fiiction force model is non-linear. 

6. The method of claim 3, wherein the slip-friction modd is defined by the following 
formulas: 



characteristics. 



4. 



The method of claim 3 wherdn the slip-fiiction model is a linear model. 




(eq. 18) 



(eq. 19) 
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where /i is a fiiction value, 2 is a slip value, is a slip slope, A„ax is a msmm^} slip 
value, Firac is a longitudinal traction force and F^is a vertical load on a wheel. 

7. The method of claim 5, wherein the non-lmear slip-fiiction model is a Continuos time 
Auto-Regressive Moving Average model (CARMA). 

8. The method of claim 5, wherein the non-linear slip-fiiction model is a peak 
parametrised model. 

9. The method of anyone of claims 3-8, further comprising a step for fiiction estimation 
(RFI), whidi provides an estimate of the maximal utili2able fiiction force Ow) and the 
initial slope (k) of the (current) slip curve, which can be used as input parameter for the 
slip-fiiction force model. 

10. The method of anyone of claims 3-9, wherein the slip referrace value is calculated 
according to the slip-fiiction force model by an adaptive filter. 

1 1 . The method of claim 10, wherein the adaptive filter is a Kahnan filter. 

12. The method of anyone of the preceding clainofi, wherein the torque control is achieved 
by fiiel blocking. 

13. The method of claim 12, characterised in that the fiiel blocking comprises a cascade 
control structure. 

14. The method of ai^one of the preceding claims, wherdn the torque control is achieved 
by controlling the throttie angle. 

15: The method of anyone of the preceding claims, wherein the torque control is achieved 
by controlling the ignition angle. 



wo 2004/067307 



-52- 



PCT/EP2003/000807 



16. The method of anyone of the preceding claims, wherein the torque control is achieved 
by controlling the wheel breaks. 

17. The method of anyone of the preceding daims, wherdn the dynamical state parameters 
comprise at least one ^eel speed parameter {o^) 

18. The method of claim 17, wherein the wheel speed parameter (a>^) is obtained by an 
anti-locking system (ABS). 

19. The method of anyone of the preceding claims, wherein the dynamical state parameters 
comprise an engine torque parameter (Te„g). 

19. The method of anyone of the preceding claims, wherein the dynamical state parameters 
comprise a vehicle velocity parameter (v^). 

20. The method of anyone of the preceding claims, wherein the dynamical state parameters 
comprise a slip curve slope parameter (k^) . 

21. The method of anyone of the precedmg claims, wherdn the slip reference vahie 
depends on the vehicle's lateral forces. 

22. A method for determinmg the traction force (Ffr«e) of a tire on a vehicle, in particular in 
combmation with one of the preceding methods, based on a slip-fiiction modd 
according to the formula 

where is a fiiction value, A is a tire slip value, Ftnci^ a traction force, Fgis the vertical 
load on a wheel and fi(X) is a slip curve according to the slip-fiiction model. 

23. The method of claim 22, wherdn the non-linear slip-fiiction modd is a Continuos time 
Auto-R^essive Moving Average model (CARMA). 
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24. The method of claim 22, wherein the non-linear slip-friction model is a peak 
parametrised model. 

25. The method of anyone of claims 22-24, wherein the slip reference value is calculated 
according to the slip-friction force model by an adaptive filter. 

25. The method of claim 24, wherein the adaptive filter is a Kahnan filter. 

26. A system for controlling the wheel slip (A) of at least one wheel on a vehicle, 
comprising: 

- a device for obtaining at least one dynamical state parameter of the vehicle, 

- a device (G(s)) for determining a slip reference value indicative of the wheel slip at 
a determined drivdine torque on the basis of a dynamic model of the relevant 
vehicle components, and 

- a device (R(s)) for controlling the driveline torque (Tj) of the vehicle by 
determining a torque manipulation signal based on the dynamical state parameters 
and the slip reference value. 

27. A computer program product including program code for carrying out a digital signal 
processing method, whai executed on a conq)uter system, for controlling the wheel slip 
(A) of at least one wheel on a vehicle, comprising: 

- a device for obtaining at least one dynamical state parameter of the vehicle, 

- a device (G(s)) for determining a slip reference value indicative of the wheel slip at 
a determined driveline tonjue on the basis of a dynamic model of the relevant 
vehicle components, and 

- a device (R(s)) for controlling the driveline torque (Z) of the vehicle by 
determining a torque manipulation signal based on the dynamical state parameters 
and the slip reference value. 
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